CoefRcient of tangential restitution for the linear dashpot model 
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The linear dashpot model for the inelastic normal force between colliding spheres leads to a 
constant coefficient of normal restitution, £„ =const., which makes this model very popular for 
the investigation of dilute and moderately dense granular systems. For two frequently used models 
for the tangential interaction force we determine the coefficient of tangential restitution et, both 
analytically and by numerical integration of Newton's equation. Although e„ =const. for the linear- 
dashpot model, we obtain pronounced and characteristic dependencies of the tangential coefficient 
on the impact velocity et = £t (g)- The results may be used for event-driven simulations of granular 
systems of frictional particles. 

PACS numbers; 45.70.-n, 45.50.Tn, 45.50.-j ,45.70.-n, 62.20.Fe 
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I. INTRODUCTION 

The interaction of dissipatively colliding spheres is de- 
scribed by the interaction force law. In the case of 
smooth spheres, that is nonfrictional spheres, the force 
is directed along the unit vector e = {ri — rj) / \ri — rj\. 
Thus, only the normal component of the impact veloc- 
ity, gn = {vi — Vj) ■ e, (with t/j/^ being the center-of-mass 
velocities of the particles i and j) changes during the con- 
tact. The normal component of the velocity after a col- 
lision is then obtained by integrating Newton's equation 
of motion for the mutual deformation ^ of the particles 



TOoffC + Fr, 



0, ^(0)=5„, e(0)=0, (1) 
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and where Fn is 



with S,{t) = max (0, Ri 
the model-specific interaction law. 

Alternatively, the collision may be described using the 
coefficient of restitution, relating the pre-coUisional and 
post-coUisional normal relative velocities. 
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with tc being the duration of the collision. Throughout 
this paper, primed variables denote post-collisional quan- 
tities. Consequently, the choice of the particular force F„ 
determines the coefficient of restitution e„. 

From its definition, Eq. obviously < e„ < 1 

where £„ = 1 describes an elastic collision. The assump- 
tion £„ =const. is widely used throughout the litera- 
ture on granular gases and other dilute granular systems. 
This assumption is not in perfect agreement with physi- 
cal reality, see [J , but it simplifies the analysis of kinetic 
and hydrodynamic equations considerably and its use is, 
therefore, justified. 

Vice versa, one may ask which force laws lead to 



En =const. For the rather general Ansatz Fn 



F^^'^ with F^f^ oc F^:^'- 



dissipative components of the interaction force, a dimen- 
sion analysis P, [3, Q shows that only combinations with 
2(c — a) +h{l + a) = lead to e„ =const. Assuming 
h ~ I, i.e., a Hnear dependence of the dissipative force on 
the velocity (which should be justified at least for small 
impact rate), we see that the choice a = 1, c = ful- 
fills the above condition. This corresponds to the linear 
dashpot model 



Fn = -kn£,~"in^- 



(3) 



Indeed, the linear dashpot model is frequently used in 
Molecular Dynamics simulations of granular systems, e.g. 
0, i, S, H, 0, H- The cocflicicnt of restitution can be 
found by integrating Eq. ([1]) with the definition Eq. ([2]), 
where the end of the collision > is given by the 
condition ^ = which takes into account that there is 
only repulsive interaction between granular particles, see 
Q for details. The coefficients of the force law, fc„ and 7„ 
translate into the impact-velocity independent coefficient 
of restitution Q 
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cx ^^"^ for the elastic and 



(5) 

There are several other force laws in the literature, 
for review see [13, [HI ■ Some of them are certainly bet- 
ter suited to describe the mechanics of colliding spheres, 
however, none of them leads to a constant coefficient of 
normal restitution. The condition £„ =const. in turn is 
essential for an entire class of scientific literature in the 
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field of dilute granular gases. Therefore, here we restrict 
ourselves to this important case. 

For the case of frictional particles, in general, a 
particle-particle interaction causes not only a change in 
the normal component of the relative velocity but also a 
change of its tangential component as well as the parti- 
cles' rotational velocity. Let us denote the relative veloc- 
ity of the particles in the point of contact by 

g = Vi-Vj- (^Ridi + i^2^^2) X e (6) 

with f2j/j and Ri/j being the angular velocities and radii 
of the two particles. Its projection to the tangential plane 
in the point of contact reads 

gt = -ex (ex g) . (7) 

Similar to the normal direction, the change of the velocity 
in tangential direction is described by the coefficient of 
tangential restitution. 



In contrast to the coefficient of normal restitution, here 
— 1 < £t < 1, that is, there are two elastic limits. The 
case e* = 1 corresponds to smooth particles, that is, the 
tangential velocity and, thus, the angular velocities of 
colliding particles do not change. The other elastic case 
et = —I corresponds to rough particles. One may think 
(in 2d) of gear wheels made of a very elastic material. 
When such particles collide, the tangential component 
of their relative velocity is reverted. The case, Et — 0, 
describes the total loss of relative tangential velocity after 
a collision. 

While the normal force is given by Eq. ^ , commonly 
in Molecular Dynamics simulations the change of the tan- 
gential velocity during an impact is described by tan- 
gential force laws which will be introduced in the next 
section. In a similar way as described above for the co- 
efficient of normal restitution, one can analyze the tan- 
gential relative velocity of colliding particles to obtain 
the coefficient of tangential restitution et- As the coeffi- 
cients of restitution are a direct consequence of the actual 
trajectory of the particles during contact, its functional 
dependence on the material properties and the impact 
velocities depends on the chosen force law. 

It is the aim of this paper to characterize the coefficient 
of tangential restitution for different expressions for the 
tangential interaction force between colliding particles as 
commonly used in Molecular Dynamics simulation. In 
particular, we are interested in the important case that 
the coefficient of normal restitution is independent of the 
impact velocity, e„ = const, as it follows from the linear 
dashpot model, Eq. 



II. TANGENTIAL FORCES 
A. Coulomb law for static friction 

The normal force between contacting spheres is deter- 
mined by their mutual compression ^ and the compres- 
sion rate ^. This is true not only for the linear dashpot 
model, Eq. ([3]), but more generally for all non-adhesive 
collision models, see [13, El: and is due to the fact that 
the interaction force is a bulk material property. In con- 
trast, the tangential force is not only a bulk property 
but also significantly determined by surface properties, 
e.g. roughness. 

The usual textbook formulation of friction distin- 
guishes between static and dynamic friction. If the par- 
ticles in contact slide on each other the friction force Ft 
(absolute value) is 

Ft = fiFN , (9) 

where Fn is the absolute value of the normal force at 
contact. Thus, the tangential force is limited by the 
Coulomb friction law, Eq. If the particles do not 

slide on each other, i.e. if the tangential relative ve- 
locity at contact is zero {gt — 0), the friction force is 
only indirectly defined. Namely, it assumes the value 
necessary to keep the particles from sliding as long as 
the resulting force does not exceed the Coulomb limit, 
Eq. ([9]). Hence, in this formulation there is no force law 
for static friction - the friction force is not determined 
from geometric properties like deformation. This makes 
the application of the Coulomb friction law in Molecular 
Dynamics simulations difficult: For the numerical inte- 
gration of Newton's equation we need in each time step 
the forces acting on the particles. These forces must be 
expressed in terms of positions and velocities of the par- 
ticles. Therefore, Coulomb's law which can a priori not 
be expressed in terms of positions and velocities must be 
modeled by a function in these variables. 

Ft = -mini ^ ^ X (10) 

[j (^Wl,W2,Wl,U2, • ■ ■ J , 

where the model specific functional /() depends on the 
history of its arguments, wi(t); < t < i, and the other 
arguments likewise. The choice of this functional is not 
unique but ambiguous to a certain degree. 

B. Model by Haff and Werner 

The most simple representation of the force scheme Eq. 
(fTU)) is the model by Haff and Werner jj] 

Ft ^ - min [fiF,„ jtgtit)] , (11) 

with the tangential component of the relative velocity at 
the point of contact, gt(i), given in Eq. (Throughout 
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this paper we call the components of the relative velocity 
before the collision g„ and gt- Time dependent veloci- 
ties that vary during the collision are called as gn (t) and 
gt{t). The final velocity components are named g'^ and 
g'f.) Without loss of generality, here and in the follow- 
ing the tangential velocity at the contact point shall be 
positive. The case of negative tangential velocity can be 
deduced by reflection and leads to identical results. 

Thus, the model by Haff and Werner assumes shear 
damping ~ gt for small velocity, limited by Coulomb's 
law. The model was successfully applied in many Molec- 
ular Dynamics simulations of granular matter although 
there appear problems when rather static systems are 
simulated, see [13] for a detailed discussion. 

However, if we consider the coefficient of normal resti- 
tution which corresponds to the model, Eq. (jlip . we 
notice more serious difficulties: both alternatively acting 
tangential forces due to Eq. pT|) cannot lead to nega- 
tive relative velocities, gt, after a collision. This is ob- 
vious for Ft oc —gt and will be shown in Sec. IIII Al for 
F oc — F„. Consequently, the model by Haff and Werner 
cannot yield coefficients of tangential restitution of neg- 
ative value. This is a serious inadequacy of this model. 
We will come back to this problem in Sec. IHI Bl 



C. Model by Cundall and Strack 

The model by Cundall and Strack [T3| mimics static 
friction by means of a spring acting in tangential direc- 
tion with respect to the contact plane. The spring is 
initialized at the time of first contact, t = 0, and ex- 
ists until the surfaces of the particles separate from one 
another after the collision. The elongation 



C{t) = J g' {t')dt' 





(12) 



characterizes the restoring tangential force, limited by 
Coulomb's law. Thus 



Ft = -min(/iF„,fctC) 



(13) 



Consequently, when fiFn < ktC, that is, the Coulomb law 
applies, the spring assumes the elongation C. = ^Fn/kt- 
The energy stored in the spring may be released in a later 
stage of the collision. From Molecular Dynamics simula- 
tion we know that this model is much better suited to de- 
scribe static behavior of granular matter As shown 
below, it yields also negative values of the coefficient of 
tangential restitution for appropriate choice of the initial 
relative velocity at the point of contact, e.g., g„ and gt- 
This is due to the fact that the internal spring acts like a 
reservoir of energy for the relative motion of the particles 
in tangential direction. In the first part of the collision 
the spring is loaded and may release the stored energy 
towards the end of the collision. This way, the tangen- 
tial component of the relative velocity may change its 



sign. The coefficient of tangential restitution which cor 
responds to the model by Cundall and Strack, Eq. 
is discussed in detail in Sec. Ill CI 



III. COEFFICIENT OF TANGENTIAL 
RESTITUTION 

A. Pure Coulomb Force 

Before discussing the most common tangential forces 
used in Molecular Dynamics simulations let us derive 
some general expressions which apply to the limiting case 
of pure sliding: Thus, for the moment we assume simpli- 
fying that during the entire collision the friction force is 
not sufficient enough to stop the tangential relative mo- 
tion. Consequently, Ft = ^Fn , from the beginning to the 
end of the collision. The results derived here are valid 
independently of the functional form of the normal force. 
Therefore, the function Fn remains unspecified, except 
for the fact that it is a function of time, defined in the 
interval (0,ic). 

During the contact the change of the velocities in nor- 
mal and tangential direction obey Newton's law. 



dt 



with 



1 

mes 
1 

m-cff 



-Fn: 



dgtit) 



dt a 



Ji 



B4 



(14) 



(15) 



which may be obtained from Eq. ([7]) with the general 
paradigm of nearly instantaneous collisions, that is, the 
unit vector e does not change during the collision. The 
validity of this approximation will be discussed briefly in 
Sec. HYl 

Using the definition of the coefficient of normal resti- 
tution, Eq. ((21), we write for the change of the normal 
component 



F.a{t)dt = (1 + £„) rricsgn ■ 



(16) 



Again we assume the tangential velocity at contact to be 
positive. During the collision the friction force assumes 
the value Ft = —fiFn{t), therefore, the tangential compo- 
nent of the relative velocity in the point of contact after 
the collision reads 



g't-gt^- r Ftdt^-^ 
a Jq a 

_ ^(H-e„)mcff 



Fndt 



(17) 



-gn 



Hence, with the definition of the coefficient of tangential 
restitution we find the relation 



et = 1 



gt 



(18) 
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always smaller than the magnitude of the shear damping 
force, 7t(?t, that is, the tangential force in the beginning 
of the collision is equal to the Coulomb force. We inte- 
grate Eq. pT|) and find that during this first stage of the 
collision, the relative tangential velocity decays as 

9t{t) = 9t-\ [cosujnt - 1) (22 

a 

If the magnitude of the shear damping force drops below 
the Coulomb force, the tangential force is governed by the 
other branch of the force law. Ft = "ftgt- The transition 
takes place at time t^^i when 



[m/s] 



(23) 



FIG. 1: The coefRcient of tangential restitution as function 
of the tangential impact velocity gt and the normal impact 
velocity <;„ in the case of pure Coulomb friction. The param- 
eters are fj, = 0.4 and kn = 10® N/m. 



independent of the functional form of the normal and 
tangential force laws (see [isl [I6;] ) . 

As an important consequence we find that the coeffi- 
cient of tangential restitution does significantly depend 
on both the normal and the tangential relative veloci- 
ties. Formally, Eq. (|18p leads to values of £4 outside 
its range of definition. However, the basic assumption of 
pure Coulomb friction implies that the particles do not 
stop sliding on each other, that is, during the entire con- 
tact the particle remain in the Coulomb regime where 
Ft — —jiFn. Hence g't>Q and, thus St > 0. 

Consequently, for pure Coulomb friction wc obtain 



St 

see Fig. [TJ 



max 0, 1 



A*(l 



I Weff gn_ 

gt 



(19) 



B. Model by HafF and Werner 

First we want to present an analytic theory for the 
elastic collision in normal direction, that is 7„ — in Eq. 

The time dependent solution of the corresponding 
Newton's equation, Eq. ([T]), reads 



^(f) = — smujnt . 



(20) 



where LUn = ykJrrhM. Note that the compression ^(i) 
is independent of the tangential relative motion of the 
particles. With Eqs. ([9]) and ([3]), the limiting Coulomb 
force is, therefore. 



IJ-Fn 



sin UJnt. 



(21) 



When the collision starts at time t = at finite tangential 
velocity 34, the magnitude of the Coulomb force fiFn is 



Inserting Eq. (1221) obtain the first switching time ts^i 
between the regimes, 



sin UJnt s,i = 7t 



aiu^id - 1) + ^a^Lul - (52 - 2Sh, 



o^'^^l + it 



6 = 



agt 



IJ^rriesgn 



(24) 



If the initial tangential velocity fulfills the inequality 



gt_ ^ fJ-mcs 



gn 



a 



a2 



1? 



(25) 



Eq. ([24)1 does not have a real solution and the Coulomb 
regime is active during the entire collision, (cf. Fig. [2l 
dashed line). In this case the coefficient of tangential 



Coulomb force 

Haff & Werner, large 

Half & Werner, medium 

Haff & Werner, [Ig^ small 




time t 

FIG. 2: Sketch of the tangential force Ft(t). When the colli- 
sion starts, the Coulomb force is always active since for finite 
gt always jtgt > jJ-Fn = 0. For large initial tangential veloc- 
ity, the condition jtgtit) < fj,F„{t) is never fulfilled during the 
entire condition. For smaller gt the tangential velocity drops 
up to t = ti,s where ytgt{tc) = l-i'Fn{tc) and the Haff-Werner 
force jtgt becomes active. The tangential velocity decays fur- 
ther and eventually at time ts,2 the Coulomb force becomes 
active again. 
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restitution reads 



2/jmoff gt 
a gn 



(26) 



which coincides with Eq. p8)) for the case e„ = 1. 

For smaller initial tangential velocities Eq. has 
a solution, that is, after the time tg^i the Haff-Werner 
branch of the tangential force becomes active and the 
tangential velocity decays exponentially: 



9t{t) =gt(ts,i)exp 



.11 
a 



(27) 



Since the tangential velocity in the Haff-Werner model 
cannot drop to zero, at some later time ts,2 > ^s,i the 
switching condition must be fulfilled again 



sin(a;„ts,2) = ltgt{tsa) 



(28) 



see Fig. [5] (dotted and dash-dotted lines) . This time is 
determined by 



exp 



It 



{tsa - ts,l) sin {uJnts,2) = gt (tss) 



^kfign 



(29) 

For t > is, 2 the Coulomb branch of the tangential force 
becomes active again for the rest of the collision. The 
final velocity reads 



g't = gt {ts,2) 



a 



cos {ujnts.2) + 1] 



and the coefficient of tangential restitution is 



£t 



El 
gt 



(30) 



(31) 



In the case of small tangential velocities, the time at the 
beginning and the end of the collision where the Coulomb 
regime is active is negligible and the decay of the tangen- 
tial velocity is mostly determined by the shear damping 
force. In this case the coefficient of restitution reads 



£t = exp 



aUJn 



(32) 



Figure [3] shows the coefficient of tangential restitution 
as function of the normal and tangential components of 
the impact velocity. We determined the coefficient of 
tangential restitution also by solving the equation of mo- 
tion, Eq. (fT4|) . numerically and analytically as described 
in this section, leading to perfect agreement. 

In the case of large tangential velocity g„, the coeffi- 
cient of restitution behaves like in the pure Coulomb case, 
that is, and for small values of gt the coefficient reaches 
the value predicted by Eq. ([32| (see Fig. [3l coarse lat- 
tices). 

In a good approximation the coefficient of tangential 
restitution in the model by Haff and Werner can conve- 
niently be described as either a non-negative constant 
or the dependence given by Eq. (llSp . whatever is larger: 



£t 



max 



2^mcff ffn 

a gt 



(33) 




[m/s] 



FIG. 3: The coefficient of tangential restitution resulting from 
the model by HafT and Werner as function of the normal and 
tangential components of the impact velocity. The coarse 
lattices display the solution for pure Coulomb force, Eq. (I19[l . 



The parameters are fc„ = 10 N/m, /i — 0.4, 74 = 10 Ns/m 
(top), 7t = 1 Ns/m (middle), 7t = 0.1 Ns/m (bottom) 



Thus, the heuristics used in [H, [T1|T3| (Eq. ((331) in a 
different notation) is justified, provided > 0. 
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C. Model by Cundall and Strack 

1. Equations of motion 

The collision model by Cundall and Strack is described 
by Eqs. d [H]) and the force laws, Eqs. The 
equation of motion reads, thus. 



(34) 



agt-Ft{C,FN) =0 
e(0) = 0; m=9n; C(0)=0; C(0) = .g* 
The tangential force 

\Ft\^ mm {ktC^f^KO (35) 
counteracts a the elongation of the tangential spring, that 



Ft = -sgnC min {ktC, M^nC) 



(36) 



The equation of motion for the tangential degree of free- 
dom thus reads 



aC + SgnC min (fctC, M^nO = • 



2. Numerical results 



(37) 



Before discussing more general properties of the col- 
lision model let us look to the typical structure of the 
coefficient of tangential restitution as it follows from the 
model by Cundall and Strack. The set of equations ([Ml 
1571) can be integrated numerically. Figure [3| shows the 
coefficient of tangential restitution as a function of the 
normal and tangential components of the impact veloc- 
ity. 

As a reference the plot on the top of Fig. [Jj shows the 
case of elastic normal force 7„ = 0. In the left panel we 
introduce damping of the normal force of different mag- 
nitude, i.e., 7„ > or £„ > 0, respectively. In the right 
panel the spring constant kt is varied. For small tangen- 
tial velocity, gt, the coefficient of tangential restitution 
approaches in all cases a constant whose value depends 
on both kt and £„. Contrary, in the case of large tangen- 
tial velocity or small normal velocity, St(gmgt) behaves 
as for the case of pure Coulomb friction. Consequently, 
for vanishing normal velocity, ffn ^ 0, see Eq. ([T5|l . 

The functional form of St{gn,gt) reveals a more com- 
plicated behavior for increasing stiffness of the tangential 
spring, kt, see Fig|4]e and f. For larger stiffness the func- 
tion St{gmgt) develops an increasing number of valleys. 
We will discuss this behavior below in Sec. IIII C 71 

In contrast, the choice of the damping in normal di- 
rection does not lead to qualitative changes of et{gmgt)- 
Only for very small coefficients of normal restitution we 
find a significant change of the form of St as compared to 
the case of pure tangential damping, e„ = 1. 




FIG. 4: The coefficient of tangential restitution as function 
of the normal and the tangential components of the impact 
velocity. The top figure shows the case of pure tangential 
damping, that is, e„ = 1. The spring constants acting in nor- 
mal and tangential direction are k„ = kt = 10®N/m. The 
left panel shows the effect of the coefficient of normal resti- 
tution on Et: a) e„ — 0.8, b) e„ — 0.6, c) £„ = 0.4. The 
right panel shows the influence of the tangential spring con- 
stant: d) kt = 0.5 A;„, e) kt — 2fc„, f) kt = 10fc„. In all plots 
the Coulomb friction parameter is ^ = 0.4. The numerical 
solution agrees perfectly with the analytical solution. 



3. Elastic normal spring 

The numerical results presented in the previous section 
suggest that the choice of the damping parameter, 7„ or 
Sn, in normal direction is of much less importance than 
the choice of the tangential spring constant fcj. It is worth 
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to discuss an approximate theory for the hmit of elastic 
interaction in normal direction, that is, 7„ = in Eq. 

For given elastic and dissipative material constants, 
kn, kt, 7n — 0, kt, and ^, the time dependence of the 
compression, ^(i), and the normal velocity, gn{t), depend 
only on the normal component of the impact velocity at 
the point of contact, = (?ri(0), but they are indepen- 
dent of the tangential component gt- Integrating Eq. 
([M)) we obtain 

^(t) ^ 5n ^^^^^ ^ ^ g ^ ^ ^ ^3g^ 

with a;„ defined in Eq. ([5]). 

Assume at time tg (Iq is not necessarily the time of im- 
pact) the tangential component of the relative velocity is 
gt (to) and the tangential motion occurs in the Cundall- 
Strack regime. At this instant the tangential velocity is 
gt{to) and the elongation of the Cundall struck spring 
C{to). The solution of the equation of motion, Eq. (HU, 
for the tangential velocity and the elongation of the tan- 
gential spring is then 



C{t) = C (to) cos [utit^to)] + 



9t (to) 



sm[LOt{t-tQ)] (39) 



gt{t) = -C (^o) sin [ujt [t - io)] + gt {to) cos [ujt {t - to)] , 

(40) 

with ujt — ^Jktfa. On the other hand, if the tangen- 
tial motion is governed by Coulomb's friction law, the 
solution of the equation of motion, Eq. p7p . is 



when the Coulomb force equals the force according to 
the Cundall-Strack force. 



[ikni(ts)^kt\^{ts)\ 



(43) 



Contrary, if at present time the motion is governed 
by the Coulomb force, determining the time when the 
Cundall-Strack regime will take over is less straight- 
foreward: If the system is in the Coulomb regime, the 
elongation of the tangential spring is determined by Eq. 
(H^ . Thus, in this regime the Coulomb force and the 
Cundall-Strack force are equal. The regime changes if 
after an infinitesimal time t -\- dt, the Cundall-Strack 
force according to Eq. p2)) exceeds fj, times the Coulomb 
force, Eq. ([55]) . that is, the next switching time is deter- 
mined by the time when the derivatives of the forces in 
both regimes equal one another. Therefore, the collision 
switches from the Coulomb regime to the Cundall-Strack 
regime at time ts with 



l^gn kn 



COs{uJnts) = gt (ts) 



(44) 



where gt{t) is governed by Eq. (j4ip . 

Finally, we have to determine whether the collision 
starts in the Coulomb regime or in the Cundall-Strack 
regime. We consider the potential change of the Cundall- 
Strack force and the Coulomb force in an infinitesimal 
time interval. The condition to start in the Coulomb 
regime reads, thus. 



fikngn < hgt 



(45) 



gt{t) = gtih) - -sgn [C {to)] / dt'F^ it') 



= gtito) + sgn[C{to)] 



a 



X [cos (cOnt) - COS (uJnto)] (41) 

and the absolute value of the spring length is 

at) = (42) 

As in the case of the Haff- Werner law, discussed in 
Sec. IIII Bl the force may change during a collision be- 
tween the Cundall-Strack regime, \Ft\ = ktC,, employing 
the elastic spring of stiffness kt and the Coulomb regime. 
Ft — iJLknS,. (For a detailed analysis of the switching 
properties see Sec. IIII C 51 ) In order to determine the 
coefficient of tangential restitution we have to determine 
the times when the regime changes and combine the cor- 
responding partial solutions, Eqs. and (|¥T|) . corre- 
spondingly. Whereas in case of the Haff- Werner force 
there are only zero or two changes of the regime, we will 
see that in case of the Cundall-Strack force there may 
occur multiple changes of the regime, see Sec. IIII C 51 

If the motion is governed by the Cundall-Strack force, 
the system changes to the Coulomb regime at time ts 



This inequality, together with Eqs. (|39II411 l43l [44)) de- 
scribe the tangential motion of the colliding spheres. We 
determined the coefficient of tangential restitution by 
combining the piecewise solutions for the Cundall-Strack 
regime and the Coulomb regime with regard to the corre- 
sponding switching times tg. The result agrees perfectly 
with the numerical solution shown in Fig. HI 

Let us now discuss the special case of pure Coulomb 
friction. If the condition if fulfilled, the dynamics 
starts in the Coulomb regime and switches to the Cundall 
Strack regime if Eq. ([44]) is fulfilled. However there may 
be no real solution of Eq. (M)) and hence the Coulomb 
regime may be active during the entire collision. In that 
case the coefficient of restitution is described by equation 
P^ . Using Eq. (|1T|) with to = and inserting this in 
Eq. (gU one finds 

— COS {uj„ts) = gt~\ cos {Ujnts) - 1] , (46) 

kt a 

assuming (without loss of generality) a positive gt yield- 
ing sgn(^ = 1. Solving this equation for cos{uJnts) one 
finds 

/ kn TOcffN / , N 9t rrics 

M l cosiunts) = M • (47) 

\kt a J gn a 
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This equation has no real solution if 
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a 





(48) 



as |cosw„ts| cannot be larger than one. Hence if both 
inequalities and ([35]) are fulfilled the dynamics stays 
in the Coulomb regime. Furthermore is possible to cal- 
culate the value of et which is reached when gt tends to 
zero. It is obvious from the inequality that the dy- 
namic starts for small gt in the Cundall-Strack regime 
and that there is a small time interval before the end 
of the collision where the dynamics is governed by the 
Coulomb regime. This time interval is proportional to 
gt/gn which is small by construction (cf. Eq. to- 
gether with Eqs. ((38l) and ((39)) ). Hence, the expression 
cos(w„to) in Eq. (gT]) is cos(w„tc) + 0{{gt/gn)^)- The 
small correction 0{{gt/ gnY) can be neglected and the fi- 
nal velocity g't can be approximated as g't = gt cos{u!ttc) 
yielding 



lim et{gn,gt) = cos(a;ttc) cos 
gt^O \ LU, 

(cf. Eq. (001) with to=0 and < = i^)- 



TTLUt 



(49) 



005 



4- Scaling properties for the case of elastic normal springs 




For the case 7„ — 0, that is, when the dissipation of 
the motion in normal direction can be neglected, £„ = 
1, apparently, the collision as described by the Cundall- 
Strack model depends on 7 parameters: m, J, fc„, kt, 
gn, and gt. By using appropriate time and length scales 
we can reduce the number of free parameters to 3. The 
length scale is the maximum compression ^max in normal 
direction: 



6r 



TOeff 



gn 



(50) 



The obvious time scale T of the problem is the duration 
of the collision, tc = tt/wq. To simplify the resulting 
expressions we drop the prefactor tt and define 



Wo 



"T-cff 
kn 



(51) 



The scaled variables are, thus, the scaled deformation, 
the scaled length of the tangential spring, and the scaled 
time, 



t/T. (52) 



Taking into account that a/nies with a defined by Eq. 
(fT5)) , for identical homogeneous spheres reduces to a pure 
number, a/uias — 2/7, the equation of motion, Eq. 
is 



i -I- a; = 

..^7 , . fkt 

z H- -sgni; mm I 7—2;, [ix 



a;(0)=0; z(0) = 0; i(0) = 1 ; i(0) 



(53) 



El 
gn 



FIG. 5: For elastic normal restitution, the Cundall-Strack 
model reveals the simple scaling form, Eq. (|53|l . Despite of 
its rather complicated functional form, the plot of et {gn,gt), 
is invariant when scaling the impact velocity. Top: velocity 
range . . . 0.02 m/sec, bottom ... 2 m/sec. The other pa- 
rameters are kn — kt — 10^ N/m, fi — 0.4. The pictures are 
identical since the coefficient of restitution depends only on 

gn/gt- 



where dots denote time derivatives with respect to the 
scaled time r. Consequently, the parameters of the sys- 
tem are /u, kt/kn, and gt/gn- As an example, the depen- 
dence on gt / gn is demonstrated in Fig. [5l When scaling 
both velocities by 100, the resulting picture is identical. 



5. Switching between friction regimes 

During a collision, depending on the impact velocity 
and the material parameters the relative motion of the 
particles at the point of contact may change its character, 
namely it may change to and fro the Coulomb regime 
where the friction force is determined (i.e. capped) by 
the normal force and the Cundall-Strack regime where 
the magnitude of the friction force is determined by the 
length of the tangential spring alone. 

This change of regime occurs also in the Haff- Werner 
model, however, there are fundamental differences: As 
shown in Sec. IIII Bl there may be only zero or two 
changes, for the Cundall-Strack model we may have mul- 
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tiple changes. This property originates from the fact that 
in the Cundall-Strack regime there is no loss of energy. 
Instead, during the Cundall-Strack regime the energy of 
the relative motion is used to load the tangential spring 
whose energy can be released subsequently, that is, the 
tangential spring acts as a reservoir of energy. The only 
way to dissipate the energy stored in the tangential spring 
is by switching into the Coulomb regime and rapidly de- 
creasing the elongation of the tangential spring due to 
decreasing normal force. In this case the energy in the 
spring cannot be fully recovered. This tangential spring 
is a great advantage of the Cundall-Strack model as its 
action may lead to a negative coefficient of tangential 
restitution which cannot be achieved by the Haff- Werner 
model. 

The number of switching events between the regimes 
as described by the criteria, Eqs. (|43)) and (|44l) and the 
initial regime given by Eq. (j45p is shown in Fig. [21 




FIG. 6: Number of switching events between the Cundall- 
Strack and the Coulomb regimes during a collision as a func- 
tion of gt/gn and kt/kn- The numbers of changes are indi- 
cated by gray shading (color online) and by the numbers in 
circles. The Coulomb friction coefficient is jj. — 0.4. The 
letters a, . . . ,o at the right side (not all are shown) refer to 
qualitatively different functions (^{t) or z{t), shown in Fig. [T] 

For any value of kt / kn there is a critical gt / gn above 
which the particle stays in the Coulomb regime (region 
"0" - black color in Fig. [5]) , corresponding to case of pure 
Coulomb friction discussed in Sec. IIII Al The boundary 
of this region is given by Eq. (|l5)) and Eq. ([^5]) 

Decreasing the tangential component of the impact ve- 
locity below the limit of pure Coulomb friction the col- 
lision will switch to the Cundall-Strack regime, at least 
for a short period of time. It may repeatedly switch back 
and forth between the two regimes, dependent on the 
impact velocity and the spring constants fc„ and kf. To 
demonstrate this effect we marked 15 points in Fig. [5] 
indicated by "a" to "o" at the right-hand side of the fig- 
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FIG. 7: Elongation of the cundall Strack Spring as funtion of 
time for kt/kn — 18 and the values gt/g-n indicated by "a" to 
"o" in Fig. [6] The Couloumb limit is shown by dashed lines. 



ure (not all are shown) which correspond to decreasing 
gt/gn for fixed kt/kn- For these parameter combinations 
we show in Fig. [7|the elongation of the tangential spring 
over time (full lines). The dashed lines show the Coulomb 
limit, i.e. iiF^/kt. Indeed, for large enough tangential 
velocity (Fig. [7^) the tangential force remains on the 
Coulomb limit for the entire duration of the collision. 

For tangential velocities slightly smaller than the limit 
given by Eq. (gS]) (region "b" in Fig. M and Fig. Uh) 
the collision switches to the Cundall-Strack regime near 
its end, that is, the friction force is sufficient to stop the 
particle once. For the chosen value of the spring con- 
stants {kt 3> kn) the period of the tangential oscillation 
is much smaller than the period of the normal oscilla- 
tion. Therefore, the tangential relative velocity of the 
particles at the point of contact is reverted. Finally, the 
Coulomb limit takes effect again, the particle stays in 
the Coulomb regime to the end of the collision. As the 
tangential motion is reverted once the coefficient of tan- 
gential restitution is negative. 

For still a little smaller velocity (region "c" in Fig. 
[S] and Fig. \7]p) the collision switches earlier from the 
Coulomb regime to the Cundall-Strack regime as its 
smaller tangential energy is dissipated earlier. Hence, 
the remaining time of contact is large enough to al- 
low not only the reversal of motion and the subsequent 
switch back to the Coulomb regime as in the case "b" 
but allows an additional switch back to the the Cundall- 
Strack regime, that is, the tangential velocity changes its 
sign back to the original direction. Finally the collision 
switches back to the Coulomb regime. Since the tan- 
gential velocity has the same sign as at the time of the 
impact, the coefficient of tangential restitution is posi- 
tive. 
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Note that the amphtude of the first tangential oscilla- 
tion (after the first switch to the Cundall-Strack regime) 
increases with decreasing tangential velocity - as the time 
of the first switch shifts towards the time of maximal 
compression. Therefore, the oscillation cannot complete 
even a half period without transiting to the Coulomb 
regime. Therefore, the number of switches is twice the 
number of zeros in the elongation of the Cundall-Strack 
spring. 

For further decreasing tangential velocity this mech- 
anism is repeated, the number of switches increases in 
steps of two. After the maximum number of switches (8) 
is achieved for the region "e" in Fig. [6] and Fig. [7^ the 
number of switches decreases. This is due to the fact that 
the amplitude of the tangential oscillation now decreases 
with decreasing tangential velocity since the first switch 
to the Cundall-Strack regime takes place before the point 
of maximal compression. The number of switches now 
decreases as the first oscillation of the Cundall-Strack 
spring may now complete more than half a period. In- 
terestingly, there is a short interval of gt where the num- 
ber of switches is back to 8 due to an additional pair of 
switches close to the end of the collision. 

For small values of the tangential velocity or small 
kt there is only one switch. The particle starts in the 
Cundall-Strack regime and performs several tangential 
oscillations without violating the Coulomb condition. 
Only at the very end of the collision the particle switches 
to the Coulomb regime. As mentioned before, this is the 
only time when energy of the tangential motion is actu- 
ally dissipated. 



6. Commensurable spring constants 

Let us discuss briefly the special case of commensurable 
tangential and normal motion that occurs when the elas- 
tic constants fe„ and kt are such that the frequency of 
the tangential Cundall-Strack spring is a multiple of fre- 
quency of the normal motion, ujt = majn, m — 1, 2, 3 . . ., 
as sketched in Fig. [H 

In this case we notice a plateau in the coefficient of 
tangential restitution as a function of the components of 
the impact velocity. Fig. [9] 

This behavior becomes clear from the sketch in Fig. 
|51 If the components of he impact velocity, gn and gt are 
such that the collision starts in the Cundall-Strack regime 
and Ut is a multiple of w„, the elongation of the tangen- 
tial spring is zero at the end of the collision. Hence, in 
contrast to non-commensurable frequencies, the system 
does not transit into the Coulomb regime close to the end 
of the collision, that is, the entire collision takes place in 
the Cundall-Strack regime. As explained above, energy is 
only dissipated in the Coulomb regime, consequently the 
coefficient of tangential restitution is £t = ±1. Whether 
the value is 1 or -1 depends only on the ratio between 
the frequencies. 
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FIG. 8: Sketch of the elongation of the tangential spring for 
commensurable frequencies, = 2tj„ and LOt = 3lj„, for the 
case that the collision starts in the Cundall-Strack regime. 



The condition for the collision to start in the Cundall- 
Strack regime is described by the inequality, Eq. ((45|) . 
therefore the boundary of the plateau is given by 



9t 



[ikn 



-9n ■ 



(55) 



7. Coefficient of tangential restitution in scaled units 



As elaborated in Sec. IIII C 4| the coefficient of tangen- 
tial restitution does not explicitly depend on the param- 
eters gt, gn, kn, and kt but only on the ratios kt / k„ and 
gt/g-n- This scaling property allows to present the coef- 
ficient of tangential restitution in a more general way it 
was shown Fig. H) Figure [TUl shows St {kt/kn,gt/gn)- 

For sufficient large tangential velocity we recover the 
pure Coulomb regime (top region in Fig. 1101 For very 
small gt/gm £t oscillates between —1 and 1. As explained 
above, in the limit of vanishing tangential velocity, gt, the 
coefficient of tangential restitution is described by Eq. 
P5|) . Expressed in terms of the ration kt/kn one obtains 



9t n 1 
^ 0, — ] ~ cos 

gn 1^1 




(56) 



et = (-1)' 



(54) 



The oscillating behavior of St {gt/ gn, kt/kn) in both di- 
rections for fixed gt/gn and varying kt/kn as well as for 
fixed kt/kn and varying gt/gn may be attributed to the 
switching between the Coulomb regime and the Cundall- 
Strack regime as discussed in Sec. IIII C 51 (see Figs. [S] 
andO. 

Note that the green line in the very left of Fig. [TOlis not 
an artifact of plotting. Here the coefficient of tangential 
restitution rises very steeply to et — > 1 since for very small 
gt during the collision whose duration is determined by 
kn the tangential spring cannot be elongated enough to 
transit into the Coulomb regime. 
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FIG. 9: The coefficient of tangential restitution as a function 
Qn and gt for commensurable frequencies of the springs in 
normal and tangential direction. In this special case, et{gn,gt) 
reveals pronounced plateaus. Top: cut = uJn, middle: Ut = 
2uj„, bottom: ujt ~ 3cj„. 



IV. VARIATION OF THE UNIT VECTOR e 
DURING THE COLLISION 
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FIG. 10: The coefficient of tangential restitution as a function 
of the ratios gt/gn and kt/kn- The Coulomb coefficient is 
^ = 0.4. 




FIG. 11: During a non-central collision the unit vector e = 
(Fi — rj) I \ri — Tj\ changes its direction. 



mi = m2 = 1 g and radii _Ri = i?2 = 4 mm. We 
computed the variation of the unit vector e during the 
collision to check whether the assumption e wconst. is 
justified. Figure [12] shows the change Ae as a function 
of the normal and tangential components of the impact 
velocity. Here Ae is defined as the angle in rad between 
the unit vector at the beginning of the collision and its 
end. 

For rather soft particles, fc„ = lO^N/m (bottom part 
of Fig fT^ the unit vector changes remarkably up to about 
30° whereas for more stiff particles, /c„ = lO^N/m (top 
part) the angle is below 1°. Hence, for sufficiently hard 
particles the assumption e=const. is justified. 



All results in the previous section were obtained 
under the assumption that the unit vector e = 
(fi — T^j) / l^i ~ '^1 keeps constant during the entire colli- 
sion. This assumption is exact only for a central collision 
when gt = 0. If the particles collide with finite tangen- 
tial relative velocity, the vector e cannot be constant but 
changes its direction, as shown in Fig. [11] 

By means of numerical integration of Newton's equa- 
tion of motion for a pair of colliding spheres of mass 



V. CONCLUSIONS 

We investigated the coefficient of tangential restitution 
for linear normal forces and two different tan gen tial force 
models - the models by Haff and Werner and by 
Cundall and Strack [l|. 

For the model by Haff and Werner, we showed that 
the coefficient of restitution is strictly non-negative. In a 
good approximation its functional form can conveniently 
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FIG. 12: Variation of the unit vector e = Ti — fj / \ fi — rj\ 
during a collision as function of gn and gt- The change Ae 
is defined as the angle in rad between the e in the beginning 
and the end of the collision. The parameters are /x = 0.4, 
mi — 7712 ~ Ig, -Ri = -R2 = 4 mm, k„ — kt = 10^ N/m (top), 
k„=kt = 10^ N/m (bottom). 



be described as either a (non-negative) constant or the 
dependence given by Eq. (fT8|) . whatever is larger (see 
Eq. ([33|l ). Thus, this model is unsuitable for describing 
collisions with negative coefficient of tangential restitu- 
tion. 

For the model by Cundall and Strack the coefficient of 
tangential restitution shows a very complex behaviour. 
For certain combinations of impact velocities and mate- 
rial parameters one may observe a negative coefficient of 
tangential restitution. By adopting suitable length- and 
time-scales one can conveniently present the velocity de- 
pendence of St by only three parameters, the friction co- 
efficient /i, the ratio of the compontents of the impact 
velocity gt / gn and the ratio of the tangential and nor- 
mal spring kt/kn, provided the dissipation of the normal 
spring can be neglected (£„ = 1). We showed that the 
latter parameter kt/kn is critical for the sign of St- 

For the limit of sufficiently large tangential velocity 
there is a universal velocity dependence of et{gn,gt) 
which is not only independent of the tangential but also 
of the normal force law. We call this limit the limit of 
pure Coulomb force. 
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